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ABSTRACT 

^ ' We present band surface photometry and major-axis kinematics of stars and 

\ ionized gas of three early-type spiral galaxies, namely NGC 772, NGC 3898 and NGC 

■"nI" ■ 7782. For each galaxy we present a self-consistent Jeans model for the stellar kine- 

matics, adopting the light distribution of bulge and disc derived by means of a two- 
dimensional parametric photometric decomposition. This allowed us to investigate the 
presence of non-circular gas motions, and derive the mass distribution of luminous and 
dark matter in these objects. 

NGC 772 and NGC 7782 have apparently normal kinematics with the ionized gas 
i-C ' tracing the gravitational equilibrium circular speed. This is not true in the innermost 

JpH| region (|r| ^ 8") of NGC 3898 where the ionized gas is rotating more slowly than the 

Q . circular velocity predicted by dynamical modelling. This phenomenon is common in 

$-H ' the bulge-dominated galaxies for which dynamical modelling enables us to make the di- 

, rect comparison between the gas velocity and the circular speed, and it poses questions 

■ about the reliability of galaxy mass distributions derived by the direct decomposition 

of the observed ionized-gas rotation curve into the contributions of luminous and dark 
matter. 

^ ! Key words: galaxies: individual: NGC 772, NGC 3898, and NGC 7782 — galax- 

' ies: kinematics and dynamics — galaxies: spiral — galaxies: formation — galaxies: 

structure 



* Based on observations carried out at the Multiple Mirror Tele- 
scope Observatory (Arizona, USA) which is a joint facility of the 
Smithsonian Institution and the University of Arizona, with the 
Isaac Newton Telescope which is operated on the island of La 
Palma by the Isaac Newton Group in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astrofi'sica de 
Canarias, and with the Vatican Advanced Technology Telescope, 
which is the Alice P. Lennon Telescope and the Thomas J. Ban- 
nan Astrophysics Facility at the Mount Graham International 
Observatory (Arizona, USA). 

t Tables 4 to 12 are only available in electronic form at the CDS 
via anonvmous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http: / / cdswcb. u-strasbg. fr / Abstract . htmj . 
■f emaihpignateiasissa. it 



1 INTRODUCTION 

In the past kinematic and dynamical studies of disc galaxies 
have been mainly focused on late-type spirals. With respect 
to lenticular and early-type spiral galaxies (defined as those 
with B jT > 0.3, following Simien & de Vaucouleurs 1986), 
the study of late-type spirals has many advantages from 
both the observational and interpretative points of view. 
Late-type spirals have small or even negligible bulges, so 
they are usually described as luminous discs embedded in 
dark matter halos. They are on average more gas-rich than 
early-type spirals. Gas is mainly confined in the galactic disc, 
considered moving at circular speed and therefore used as 
direct tracer of the gravitational potential. Since the line-of- 
sight velocity distribution (LOSVD hereafter) of ionized or 
neutral hydrogen is easier to measure than the stellar one, 
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the number of spirals of known gaseous kinematics increased 
faster than that of the galaxies for which stellar motions 
were measured. These facts made late-type spirals easy tar- 
gets, particularly in the search for dark matter (see Sofue & 
Rubin 2001 for a review). 

Only recently Heraudeau & Simien (1998) and Her- 
audeau et al. (1999) measured the stellar rotation curves 
and velocity-dispersion profiles along the major axis of a 
large sample of spirals (more than 60 objects) ranging from 
Sa to Sc. In last years systematic surveys have been started 
to derive the stellar and gaseous kinematics in SO's (Bertola 
et al. 1995; Fisher 1997; Bettoni & Galletta 1997) and in 
early-to-intermediate type spirals (Corsini et al. 1999; Vega 
Beltran 1999; Vega Beltran et al. 2000). 

If the direct comparison of stellar and gaseous kinemat- 
ics allowed to unveil also in disc galaxies the presence of 
kinematically-decoupled components (see Bertola & Corsini 
1999 for a recent review), the application of dynamical mod- 
els showed that in the bulge-dominated region the ionized- 
gas velocity may fall below the circular speed (Fillmore, 
Boroson & Dressier 1986; Kent 1988; Kormendy & Westp- 
fahl 1989). These 'slowly-rising' rotation curves have been 
interpreted as the signature of the presence of pressure- 
supported gas (Bertola et al. 1995; Cinzano et al. 1999), 
which has been recently expelled by stars of the bulge but 
not yet heated to the virial temperature of the galaxy. This 
discovery poses new questions about the reliability of galaxy 
meiss distributions derived by the decomposition of the ob- 
served ionized-gas rotation curve into the contributions of 
luminous and dark matter. In fact the inner gradient of the 
gas velocity curve is usually used to fix the amount of lu- 
minous matter, which results to be underestimated if gas 
velocity rises slowly than circular speed. 

In this paper we present a study of the ionized-gas and 
stellar kinematics for 3 Sa - Sb galaxies, namely NGC 772, 
NGC 3898, and NGC 7782. We apply a self-consistent dy- 
namical model based on a two-dimensional photometric de- 
composition method to explain the different kinematical be- 
haviour of gas and stars. The goal of this work is to derive 
the mass distribution of these objects and investigate the 
possible presence of non-circular gas motions in the bulge 
region. 

The paper is organized as it follows. In Sect. ^ we give 
an overview of the properties of NGC 772, NGC 3898 and 
NGC 7782; in Sect. ^ we present the broad and narrow-band 
imaging and major-axis kinematics of stars and ionized gas; 
in Sect. ^ we discuss the photometric and dynamical tech- 
niques adopted to investigate the light and mass distribution 
of these galaxies; in Sect, fa we apply these techniques to the 
galaxies, deriving their structural and kinematic parameters; 
finally, Sect. ^ is devoted to the discussion of the results and 
conclusions. 



2 GLOBAL PROPERTIES OF NGC 772, NGC 
3898 AND NGC 7782 

The galaxies studied in this paper are a subset of the 20 disc 
galaxies observed by Vega Beltran et al. (2000). All these 
galaxies are bright [Bt < 13.5) and nearby objects (V© < 
5800 km s~^) with an intermediate-to-high inclination (i > 



45°) and their Hubble morphological types run from SO to 
Sc. 

Since the presence of ionized-gas supported by non- 
circular motions have been observed only in lenticular or 
bulge-dominated spiral galaxies, we focused our attention 
on SO - Sb galaxies. We selected unbarred galaxies accord- 
ing to the classification of both Sandage & Tammann (1981, 
hereafter RSA) and de Vaucouleurs et al. (1991, hereafter 
RC3). To better disentangle the contributions of the bulge, 
the disc and eventually of the dark matter halo to the total 
mass of the galaxy we choosed between the sample objects 
those with the more extended ionized-gas and stellar kine- 
matics (in units of i?25). At the end of this selection process 
we remained with five galaxies: an SO (NGC 980), two Sa 
(NGC 772 and NGC 5064) and two Sb spirals (NGC 3898 
and NGC 7782). 

All the selected galaxies display a smooth and regular 
morphology, except for NGC 772 which is characterized by 
a strong lopsidedness (Fig. However the symmetry of ve- 
locity curves and velocity-dispersion profiles of both ionized 
gas and stars (Fig. ^) suggests us that also NGC 772 has 
an axisymmetric structure at least in the radial region were 
the kinematic parameters were measured (corresponding to 

0.2ii25). 

An accurate analysis of the photometric and kinematic 
properties of the selected galaxies (see Vega Beltran 1999 
for details) showed that the surface brightness of NGC 980 
is characterized by a strong twisting of the isophotes and 
that the LOSVD profiles of NGC 5064 have a strongly non- 
Gaussian shape. We interpreted these two phenomena as due 
to the misalignment between bulge and disc of NGC 980 (or 
even to the presence of a triaxial bulge) and to the possi- 
ble coexistence in NGC 5064 of two counterrotating com- 
ponents, respectively. Therefore the two galaxies were dis- 
carded and are not studied here. 

An overview of the properties of the remaining galaxies, 
namely NGC 772, NGC 3898 and NGC 7782 which will be 
the subject of our investigation, is given in Table ^. Their 
available photometric and kinematical data are listed here 
briefiy. 



2.1 NGC 772 

Surface photometry of NGC 772 was obtained in the B band 
by Lu (1998), who also derived the photometric parameters 
of the exponential disc, and in the V and / bands by Her- 
audeau & Simien (1996). NGC 772 belongs to the sample of 
34 early-to-late spiral galaxies, whose major-axis stellar ve- 
locities has been recently measured by Heraudeau & Simien 
(1998). NGC 772 was observed in the 21-cm line of neutral 
hydrogen by Rhee & van Albada (1996), who obtained the 
Hi position-velocity map, global velocity profile and radial 
surface density distribution. These data show the asymmet- 
ric distribution of the Hi at radii larger than 5' from the 
centre. On the NW side a low rotation-velocity [Vhi — 100 
km s~^) component can be traced out to 9', while on the 
SE the H I emission extends to 5' with a rotation velocity of 
260 kms"\ 
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Table 1. Optical and radio properties of NGC 772, NGC 3898 and NGC 7782 



Parameter 


JNGC 772 


NGC 3o9o 


JNGC 7782 


Other name 


UGC 1466; PGC 7525 


UGC 6787; PGC 36921 


UGC 12834; PGC 72788 


Morphological type 


Sb(rs)I Sb SAS3 


Sal "'; Sa SAS2 


Sb(s)I-n ^; Sb *>; SAS3 


TT 1* J. ■ J. • 1 "J. /"l — 1 \ c\ 

Heliocentric systemic velocity (km s j 


2470 ± 10 


1184 ± 10 


5430 ± 10 


Distance (Mpc) 


34.7 


17.1 


75.2 


Major-axis position angle ° 


130° 


107° 


175° 


Apparent isophotal diameters ° 


7^2 X 4.'3 


4.'4 X 2^6 


2'.'4 X 1^'3 


Inclination 


54° 


54° 


58° 


Apparent Vt magnitude (mag) ° 


10.31 


10.70 


12.23 


Total {B — V)t color index (mag) ° 


0.78 


0.90 


0.85 


Total corrected V luminosity L,,o (Lq v)'^ 


11.2 • 10i° 


1.6 • 10i° 


9.7-1010 


Hi lincwidth at 20% of the peak (km s ^) 


473 « 


504 f 


570 s 


Hi lincwidth at 50% of the peak (km s"^) 


410 " 


469 ' 


559 8 


Mass of neutral hydrogen Mm (Mq) 


25.6 ■ 10^ « 


2.6 • 10-' f 


14.5 ■ 10^ s 


Mass of cool dust Mj (Mq) ^ 


2.3 • 10^ 


0.07- 10'^ 


2.3 • lO'' 



from RSA. 

from Nilson (1973, hereafter UGC). 

from RC3. The apparent isophotal diameters are measured at a surface brightness level of = 25 mag aresec^^. 

from this paper. The distance is derived as Vo/Hq with Vb the velocity relative to the centroid of the Local Group 
obtained from the heliocentric systemic velocity as in RSA and Hq = 75 km s^^ Mpc~l. The inclination i is derived as 
cos^ i = (q^ — gQ)/(l — Qq), where the observed axial ratio is taken from RC3 and an intrinsic flattening of go = 0.11 has 
been assumed following Guthrie (1992). 

" from Rhee & van Albada (1996), the neutral hydrogen mass has been scaled for the adopted distance. 

^ from van Driel & van Woerden (1994), the neutral hydrogen mass has been scaled for the adopted distance. 

S from Krumm & Salpeter (1980), the neutral hydrogen mass has been scaled for the adopted distance. 

^ derived following Young et al. (1989) from the IRAS flux densities at 60 and 100 urn (Moshir et al. 1990). 



2.2 NGC 3898 

Surface photometry of NGC 3898 is available in B (Barbon, 
Benacchio & Capaccioli 1978; Boroson 1981), V (Watan- 
abe 1983; Kodaira, Okamura & Ichikawa 1990; Heraudeau & 
Simien 1996), I (Heraudeau & Simien 1996), r (Kent 1988), 
J (Giovanardi & Hunt 1996; Moriondo, Giovanardi & Hunt 
1998a), H (Giovanardi & Hunt 1996), and K band (Gio- 
vanardi & Hunt 1996; Moriondo et al. 1998a). Whitmore, 
Rubin & Ford (1984), Fillmore et al. (1986), Heraudeau et 
al. (1999) measured the major-axis stellar velocity curve and 
velocity-dispersion profile of NGC 3898. The ionizcd-gas ro- 
tation curve was obtained along the galaxy major axis by 
Rubin et al. (1985) and Fillmore et al. (1986). The distribu- 
tion and velocity field of Hi were studied in detail by van 
Driel & van Woerden (1994), who also derived using the 
maximum-disc hypothesis (van Albada & Sancisi 1986) the 
mass contribution of bulge, disc and dark halo by fitting the 
combined Ha/Hi rotation curve using the photometric pa- 
rameters of bulge and disc obtained from Watanabe's (1983) 
surface-brightness profile. Other mass models for NGC 3898 
have been obtained by Fillmore et al. (1986) using both gas 
and stellar kinematics and by Kent (1988) and Moriondo, 
Giovanardi & Hunt (1998b) who adopted the ionized-gas 
kinematics by Rubin et al. (1985). 



2.3 NGC 7782 

Surface photometry of NGC 7782 has been obtained in the 
V band (Kodaira et al. 1990), the r band (Courteau 1996) 
and the H band (Moriondo et al. 1999). The only available 



bulge-disc decomposition for NGC 7782 is that of Baggett, 
Baggett & Anderson (1998) based on the data of Kodaira et 
al. (1990). No spatially resolved kinematics for the gaseous 
and stellar components have been obtained for this galaxy 
either at optical or at radio wavelengths. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Long-slit spectroscopy 

The spectroscopic observations of NGC 772, NGC 3898 and 
NGC 7782 were carried out in two difi^erent runs during Oc- 
tober and December 1990 at the 4.5-m Multiple Mirror Tele- 
scope (MMT) on Mt. Hopkins (Arizona, USA). The 1200 
grooves mm~^ grating blazed at 5767 A was used in the 
first order in combination with a l'.'25 x 3!0 slit and the Lo- 
ral 1200 X 800 CCD with pixels of 15 x 15 ^m^. It yielded a 
wavelength coverage of 650 A between 4850 and 5500 A with 
a reciprocal dispersion of 54.7 A mm~^. No on-chip binning 
was performed and every spectrum pixel corresponded to 
0.82 A by 0'.'30. 

NGC 3898 was also observed at the Isaac Newton Tele- 
scope (INT) in La Palma (Spain) on March 19, 1996. The 
Intermediate Dispersion Spectrograph (IDS) was used with 
a l'.'9 X 4!0 slit, the 500 mm camera, the AgRcd collimator, 
and tlie H1800V grating with 1800 grooves unn^^ at first 
order. This instrumental set-up yielded a wavelength cover- 
age of 240 A between 6650 A and 6890 A with a reciprocal 
dispersion of 9.92 A mm~^. No on-chip binning was applied 
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Table 2. Log of the spectroscopic observations 



Object 


Date 


Telescope 


¥\ 


P.A. 

ri 


NGC 772 


22 Oct 1990 


MMT 


3600 


130 


NGC 3898 


18 Dec 1990 


MMT 


3600 


107 




19 Mar 1996 


INT 


4 X 3600 


107 


NGC 7782 


22 Oct 1990 


MMT 


3600 


175 



on the adopted 1024x1024 TK1024A CCD. Each 24 x 
24 /xm spectrum pixel corresponds to 0.24 A by 0'.'33. 

At the beginning of each exposure, the sht was centred 
on the galaxy nucleus using the guiding TV camera and 
aligned along the galaxy major axis. The details on the slit 
position and spectra exposure times are given in Table ^. At 
the MMT spectra of a number of late-G and early-K giant 
stars were obtained with the same set up to serve as tem- 
plates in measuring the stellar kinematics. In all observing 
runs, comparison exposures of the arc lamp were taken be- 
fore and after each object integration to allow an accurate 
wavelength calibration. Quartz-lamp and twilight-sky flat 
fields were used to map pixel-to-pixel sensitivity variations 
and large-scale illumination patterns. The seeing during the 
observations was typically between l'.'2 and I'.'S FWHM. 

Using standard MIDAsQ routines the spectra were bias 
subtracted, flat-field corrected and wavelength calibrated. 
Cosmic rays were identified by comparing the counts in 
each pixel with the local mean and standard deviation (as 
obtained by the Poisson statistics of the photons know- 
ing the gain and readout noise of the detector), and cor- 
rected by interpolating. The instrumental resolution was 
derived by measuring the Gaussian FWHM of a dozen of 
unblended arc-lamp lines distributed over the whole spec- 
tral range of a wavelength-calibrated comparison spectrum. 
We found a mean value of FWHM = 2.24 ± 0.26 A and of 
FWHM = 2.57 ± 0.11 A for the MMT spectra obtained in 
October 1990 and December 1990, respectively. They corre- 
spond to an instrumental velocity dispersion of cjinatr = 55 
km s"^ and (Jinstr = 64 km s"^ at 5150 A. For the INT 
spectra we measured a mean FWHM = 0.87 ± 0.04 A that, 
in the range of the observed gas emission lines, corresponds 
to an instrumental velocity dispersion of Oinstr = 17 km s~^. 

The stellar kinematic parameters were measured from 
the absorption lines present on MMT spectra using the 
Fourier Correlation Quotient method (Bender 1990), as ap- 
plied by Bender, Saglia & Gerhard (1994). The spectra of 
the G5III star HR 7778 and K2III star HR 6415 provided 
the best match to galaxy spectra obtained in October and 
December 1990, respectively, so they were used as templates 
to measure the galaxy stellar velocities in the two runs. The 
stellar kinematics of N GC 7 7 2, NG C 38 98 an d NGC 7782 
are discussed in Sects. ^.l.l| , |5.2.l| , and ^.3.l| , and the key 
parameters are tabulated in Tabs. 4, 5 and 6. Each table 
provides the radial distance from the galaxy centre in arc- 



sec, the observed heliocentric velocity and the velocity dis- 
persion in km s^^, and the Gauss-Hermite coefficients /13 
and hi. 

The ionized-gas kinematics was derived by measuring 
the position and width of the [O iii] A5006.8 and Ha emis- 
sion lines in the MMT and INT spectra, respectively. Using 
the MIDAS package ALICE we fitted interactively a Gaus- 
sian to the emission line and a polynomial to its surrounding 
continuum. The central wavelength and FWHM (corrected 
for instrumental FWHM) of the fitting Gaussian were con- 
verted into radial velocity and velocity dispersion, respec- 
tively. The resulting velocities were corrected to the helio- 
centric frame of reference. At radii where the intensity of 
the relevant emission was low, we averaged from 3 to 7 ad- 
jacent spectral rows to improve the signal-to-noise ratio of 
the line. The ionized-gas kinematics o f NGC 772 , NGC 3898 
and NGC 7782 are described in Sects 



5.1.1 



5.2.1 



3 MIDAS is developed and maintained by the European Southern 
Observatory 



and 5.3 

and the key parameters are tabulated in Tabs. 7, 8 and 9. 
Each table provides the radial distance from the galaxy cen- 
tre in arcsec, the [O iii] A5006.8 (and the Ha only for NGC 
3898) observed heliocentric velocity and velocity dispersion 
in km s~^. 



3.2 Broad-band imaging 

The broad-band imaging of the three galaxies was carried 
out at the 1.83-m Vatican Advanced Technology Telescope 
(VATT) operated in the Mt. Graham International Obser- 
vatory (Arizona, USA) in two observing runs on March and 
November 1997. A back-illuminated 2048 x 2048 Loral CCD 
with 15 X 15 ^ra? pixels was used as detector at the aplanatic 
Gregorian focus (f/9) of the telescope. It yielded a field of 
view of 6^4 x 6!4 with an image scale of 0'.'4 pixel" ^ after a 
2x2 pixel binning. The gain and the readout noise were 1.4 
e~ ADU"'^ and 6.5 e~, respectively. 

At regular intervals during each night, different bias 
frames (typically ten) were taken to check possible slight 
bias-level variations. All galaxies were observed twice in the 
Johnson V band. The date and the duration of all the ex- 
posures are given in Table |^. A number of twilight sky fiats 
were taken at the beginning and at the end of the nights. 
No photometric standard was observed. 

The data reduction was carried out using standard 
iraf|3 routines. All the frames were bias subtracted and 
corrected for pixel-to-pixel intensity variations by using a 
mean flat fleld for each night. The different frames of each 
galaxy were shifted and aligned to an accuracy of a few 
hundredths of a pixel using common field stars as refer- 
ence. After checking that their point spread functions (PSF) 
were comparable, the frames were averaged to obtain a sin- 
gle y— band image. The cosmic rays were identified and re- 
moved during the averaging routine. Two-dimensional Gaus- 
sian fits to the field stars in the resulting images yielded the 
final FWHM measurement of seeing PSF listed in Table || 
In each final frame the mean value of the sky level was de- 
termined in a large number of 5 x 5 pixel areas. These areas 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories which are operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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Table 3. Log of the broad and narrow-band observations 



Object 


Date 


Filter 


texp 
l°l 


Seeing^ 
L 1 


NGC 772 


02 Nov 1997 


V 


2 X 180 


1.7 




07 Nov 1997 


r6450 


4 X 600 


1.9 




07 Nov 1997 


r6630 


4 X 600 


1.8 


NGC 3898 


11 Mar 1997 


V 


2 X 60 


1.4 




11 Mar 1997 


r6450 


4 X 540, 2 X 600 


1.3 




11 Mar 1997 


r6580 


2 X 500, 2 X 540 


1.3 


NGC 7782 


02 Nov 1997 


V 


2 X 360 


2.2 




10 Nov 1997 


r6450 


3 X 600 


1.8 




10 Nov 1997 


r6680 


3 X 600 


2.2 



5.2.2| , and p.3.2| , and tabulated in Tabs. 10, 11, and 12. Each 
table provides the isophotal semi-major axis in arcsec, sur- 
face brightness in mag arcsec"^, ellipticity, position angle in 
degrees, and the cos 46 Fourier coefficient. 



3.3 Narrow-band imaging 

The narrow-band Ha imaging of NGC 772, NGC 3898 and 
NGC 7782 was performed at the VATT in the same observ- 
ing runs, during which we obtained the V— band images. 

For each galaxy we obtained at least two emission- 
band and two continuum-band images using the set of in- 
terference filters kindly provided by R.C. Kennicutt. The 



10 



11 



12 



a 1 J 



> 



Seeing FWHM measured on the final resulting frame. 



were selected in empty regions of the frames, which were free 
of objects and far from the galaxy to avoid the contamina- 
tion of the light of field stars and galaxies as well as of the 
target galaxy itself. The sky value of each frame is the av- 
erage of these mean values. For the estimate of error in the 
sky determination we adopted half the difference between 
the maximum and minimum of the mean values obtained 
for the small areas. 

For each galaxy we derived a 'luminosity growth curve' 
by measuring the integrated magnitudes within circular 
apertures of increasing radius by means of the IRAF task 
ELLIPSE within the STSDAS package. Absolute calibration 
was performed by fitting the constant portion of the growth 
curves (Fig. ^) to the corresponding total magnitude Vr 
given by RC3. 

Before fitting ellipses to the galaxy images, we re- 
moved the field stars using the IRAF routines within the 
DAOPHOT package. Special care was taken to remove sat- 
ured stars and field galaxies which were edited out by hand, 
and removed by replacing them with the local average of 
counts. When bright stars were too close to the galaxy their 
light was masked out to where the galaxy light started to be 
dominant. For each galaxy we derived the surface-brightness 
profile and isophotal shape parameters by ellipse fitting to 
its isophotes using the isophote fitting program ELLIPSE 
(see Jedrzejewky 1987 for details on the fitting procedure). 
The resulting surface brightness, ellipticity, position angle 
and cos 46 Fourier coefficient radial profiles for NG C 772 , 
NGC 389 8, an d NGC 7782 are presented in Sects. ^.L2| , 

10 



14 

15 
16 
17 




50 



1 00 1 50 

r [arcsecj 



200 



Figure 1. Calibrated luminosity growth curves of NGC 772 
(filled squares, Vr = 10.31 mag), NGC 3898 [open diamonds, 
Vt = 10.70 mag) and NGC 7782 [open circles, Vt = 12.23 mag). 
For each galaxy the dashed line shows the asymptotic magnitude 
taken from RC3. 



emission-band images were taken with three different inter- 
ference filters (Ac = 6630 A, Ac = 6580 A, and A^ = 6680 
A respectively for NGC 772, NGC 3898 and NGC 7782; 
AApwHM = 70 A) in order to isolate the spectral region 
characterized by the redshifted Ho? and [N ii] (AA6548.0, 
6583.4 A) emission lines according to galaxy systemic ve- 
locity. The continuum-band images were taken through an 
interference filter (Ac = 6450 A; AAfwhm = 70 A), which 
has been selected to observe an emission-free spectral region 
sufficiently near to that of the emission-band filters in order 
to subtract off the stellar continuum in the emission-band 
images. Duration, date and filter of all the in-band and off- 
band exposures are collected in Table ^ Flat-field exposures 
of the twilight sky were taken at dawn and sunset for each 
fiher. 

The narrow-band images of the target galaxies were 
bias subtracted, flat-field corrected, aligned, averaged and 
cleaned of cosmic ray events in the same way as for the 
broad-band images. As result of this standard data reduction 
we obtained for each galaxy a single emission-band image 
and a single continuum-band image. The seeing FWHM's 
of these images as measured by fitting a two-dimensional 
Gaussian to the field stars are listed in Table ^ For each 
galaxy the best-seeing image was convolved with a Gaussian 
PSF to yield the same PSF FWHM of the worst-seeing im- 
age. Finally the continuum-band image was suitably scaled 
and subtracted from the emission-band image to obtain a 
continuum- free map of the Ha-|-[N ii] emission of the galaxy. 
The scale factor for the continuum-band image was esti- 
mated by comparing the intensity a large number of field 
stars in the two bandpasses. The resulting continuum-band 
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and continuum- free images o f NGC 772 , N GC 38 98 and 
NGC 7782 are given in Sects. ^.1.3[ |5.2.4 and ^.3.3 



4 STRUCTURE MODELS 

In order to describe the mass structure of each galaxy, we 
have to take advantage of both the photometric and kine- 
matic data available. 

We describe each galaxy by the superposition of two 
different components, namely a bulge and a disc. Both the 
components are assumed to be oblate, with the isodensity 
surfaces being similar concentric spheroids. In this frame- 
work, what we call a 'disc' is in fact a spheroidal component 
with a very high ellipticity, which we assume to be constant 
with radius. 

Once the luminosity density parameters are constrained 
by the photometry (Sect. ^^), we use the kinematics to eval- 
uate the mass-to-light ratios and to discuss the gas velocity 
distribution. 



4.1 Photometric decomposition 

Conventional bulge-disc photometric decompositions based 
on elliptically averaged surface-brightness profiles are sub- 
ject to strong systematic errors due to the different intrinsic 
shapes of bulge and disc and to the viewing angle of the 
galaxy (see Byun & Freeman 1995 for an extensive discus- 
sion). For this reason we applied to the V'— band images 
of NGC 772, NGC 3898 and NGC 7782 a two-dimensional 
parametric decomposition technique. The bulge and disc 
photometric parameters used as the initial trial in modelling 
the observed kinematics were derived by fitting directly the 
galaxy surface-brightness distributions. Our decomposition 
method is based on the technique developed by Byun & 
Freeman (1995) which we improved by introducing convo- 
lution with the seeing PSF and by weighting the observed 
surface brightness measured in each pixel according to the 
associated Poissonian noise and CCD readout noise. 



4-1.1 Two-dimensional bulge-disc parametric 
decomposition 

We assumed the galaxy surface-brightness distribution to be 
the sum of the contributions of an oblate bulge and an in- 
finitely thin disc. We adopted the r^'^* law (de Vaucouleurs 
1948) to describe the surface brightness of the bulge com- 
ponent 



Ibulge {X, y) 

h exp < 



-7.67 



' ^x^+y-^{b/a) 



-2 

bulge 



1/4 



(1) 



where le and r^ are respectively the effective surface bright- 
ness and effective radius of the bulge, x and y are the ap- 
parent distances from galactic centre along the major and 
minor axis, and {h/a)huigi: is the apparent axial ratio of bulge. 
Bulge isophotes are ellipses with constant apparent elliptic- 
ity ebulge = 1 — {b /a) bulge- 



The disc component is assumed to have an exponential 
surface brightness (Freeman 1970) 



Imsc {x, y) = lo exp 



Vx^ + y^{b/ci)aL 



(2) 



where and h are respectively the disc central surface 
brightness and scale length, x and y the apparent distances 
from galactic centre along the major and minor axis, and 
{b/ a) disc is the apparent axial ratio of disc. Disc isophotes are 
ellipses with constant apparent ellipticity edisc = l — {b/a)disc 
and the disc inclination is i = arccos {b/a)disc- 

To derive the photometric parameters of the bulge (7e, 
re and {b/a)i,uige) and the disc (7o, h and {b/a)disc) we fit- 
ted iteratively a model of the surface brightness to the ob- 
servations using a non-linear minimization based on the 
Levenberg-Marquardt method (e.g. Bevington & Robinson 
1992; Press et al. 1996). 

The galaxy centre, the bulge and disc position angles, 
as well as the sky level are not free parameters in our de- 
composition. The galaxy centre is assumed to be the mean 
centre of the ellipses which fit the galaxy isophotes. Bulge 
and disc are assumed to have the same position angle which 
is the mean value of the major-axis position angle of ellipses 
which fit the galaxy isophotes. The sky level is estimated 
as described in Sect. |^. The seeing effects were taken into 
account by convolving the model image with a circular Gaus- 
sian PSF with the FWHM measured using the stars in the 
image (Table |^. The convolution was performed as a prod- 
uct in Fourier domain before the minimization. 

For each pixel {x, y), the observed galaxy photon counts 
Tgai are compared with those predicted from the model 
Ibuige+ldise- Each pixel is weighted according to the variance 
of its total observed photon counts due to the contribution 
of both galaxy and sky, and determined assuming photon 
noise limitation. Therefore we can write 



N M 2 

[Jmod [x,y) -Igai [x,y)\ 

X — l y—1 



Igai {x, y) + Isky {x, y) 



(3) 



with X and y ranging over the whole A'^ x M pixel image. 

To derive the six free parameters of the model surface- 
brightness distribution, we adopted as the initial trial for 
minimization the values obtained by performing a stan- 
dard photometric decomposition with a parametric tech- 
nique similar to that adopted by Kent (1985). In fact we de- 
composed the observed surface-brightness profile along both 
the major axis (obtained by fitting ellipses to isophotes) and 
the minor axis (obtained from the major-axis profile scaled 
by the factor 1 — e) as the sum of an r^/^ oblate bulge plus an 
exponential infinitely thin disc. To take into account the see- 
ing effect, we truncated the major and minor-axis profiles at 
a radius corresponding to a couple of PSF FWHM from the 
centre. We assumed the minor-axis profiles of bulge and disc 
to be the same as the major-axis profile scaled respectively 
by a factor 1 — etuige = {b/a)tuige and 1 - edisc = {b/ a) disc- A 
least-squares fit of the model to the photometric data pro- 
vided the values of 7e, r^, {b/a)buige, lo, h and {b/a)dise to be 
used as the initial trial parameters for the two-dimensional 
photometric decomposition. 

The photometric parameters of the bulge and the disc 
obtained for NGC 772, NGC 3898 and NGC 7782 are given 
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Table 13. Parameters from the bulge-disc photometric decomposition 



Object Bulge Disc Fit^ 

[mag arcsec"^] ["] [kpc] [mag arcsec"^] ["] [kpc] [°] 



NGC 772 


22.7 


60.5 


10.2 


0.66 


21.5 


96.3 


16.2 


0.39 


67 


ID 




22.8 


65.4 


11.0 


0.77 


21.3 


67.0 


11.3 


0.62 


52 


2D 


NGC 3898 


20.9 


25.6 


2.1 


0.64 


21.6 


45.9 


3.8 


0.37 


68 


ID 




20.6 


18.9 


1.6 


0.71 


20.4 


29.0 


2.4 


0.52 


59 


2D 


NGC 7782*^ 


20.2 


3.6 


1.3 


0.85 


20.5 


23.3 


8.5 


0.41 


66 


ID 




19.9 


2.9 


1.1 


0.71 


20.3 


21.6 


7.8 


0.56 


56 


2D 



ID = parametric fit on the ellipse-averaged profiles, 2D = two-dimensional parametric fit. 
^ For NGC 7782 the best-fit to the observed surface-brightness distribution has been obtained adopting an exponential 
bulge (see discussion in Sec. 5.3.4). 



in Table The model surface brightnesses obtained for the 



sample galaxies are discussed in Sects. ^t!2 , ^.2.2 



and 5.3.2 



4-1.2 Test with model galaxies 

To test the reliability and accuracy of our two-dimensional 
technique for bulge-disc photometric decomposition we ap- 
plied the decomposition program to a set of artificial disc 
galaxies. We generated 100 images of galaxies with an r^''* 
oblate bulge plus an exponential infinitely thin disc. The 
scale surface-brightness, scale length, and apparent axial ra- 
tios of bulge and disc of the artificial galaxies were randomly 
chosen in the range of values observed by Kent (1985) for a 
large sample of SO and spiral galaxies. The adopted ranges 
are 



18 < /ie < 24 mag arcsec , 
3 < r-e < 10 kpc, 

0.6 < {b/a)l,nlge < 1, 

for the bulge parameters, and 

19 < /ie < 23 mag arcsec"'^, 
3 < r-e < 15 kpc, 

0.2 < (b/a)*.c < 1 



(4) 
(5) 
(6) 

(7) 
(8) 
(9) 



for the disc parameters. The parameters of the artificial 
galaxies have also to satisfy the following conditions 



{h/ a) disc < {b/ a) bulge 

< B/T < 0.8. 



(10) 
(11) 



All the simulated galaxies were assumed to be at the 
Virgo Cluster distance (d = 17 Mpc, Freedman et al. 1994) 
corresponding to a scale of 82.4 pc arcsec"'^. The pixel scale 
used was 0'.'4 pixel"^ and the CCD gain and readout noise 
were respectively 1.4 e~ ADU"'^ and 6.5 e~ in order to sim- 
ulate the VATT observations. We fixed the seeing FWHM at 
2'.'0. We added a sky background level of 10 counts pixel"^ 
which corresponds to a sky surface brightness of 21.8 mag 
arcsec"^ in the 1/— band adopting the absolute calibration 
derived from our real observations. This is a typical surface 
brightness value for a dark sky as reported by Binney & 



Merrifield (1998). An appropriate level of noise was added 
to the resulting image to yield a signal-to-noise ratio similar 
to that of the photometric data we obtained for NGC 772, 
NGC 3898 and NGC 7782 during the VATT observing runs. 

The images of the artificial galaxies have been analyzed 
as if they were real. The two-dimensional parametric decom- 
position has been applied using as initial trial parameters the 
le, Te, {h/a)buige, h, h, aud {b/a)disc values obtained with a 
standard parametric decomposition of the major and minor- 
axis surface brightness profiles measured by fitting ellipses 
to the galaxy isophotes. The fitting algorithm was able to re- 
cover the input parameters within the estimated errors with 
a scatter consistent with the results of the similar double- 
blind tests performed by Schombert & Bothun (1987), Byun 
& Freeman (1995) and Wadadekar, Robbason & Kembhavi 
(1999). 

4.2 Dynamical models 

In order to investigate the gas and stellar kinematics we use 
the self-consistent dynamical models by Pignatelli & Gal- 
letta (1999). We give here just a brief summary of the tech- 
nique itself and of the general assumptions used. 

The galaxy can be described by superposition of differ- 
ent components. For each component, we separately assume: 

- the density distribution is oblate, without triaxial struc- 
tures; 

- the isodensity surfaces are similar concentric spheroids; 

- the surface density profile follows a simple r^^'^ or an 
exponential law; 

- the velocity distribution is locally Gaussian; 

- the velocity dispersion is isotropic (ct,. = (tb = o^); 

- the mass-to-light ratio is constant with radius. 

Our model does not consider the possible presence of tri- 
axial structures (bar; triaxial bulge; tilted component; warp) 
or of anisotropy in the velocity distribution (with predomi- 
nance of radial or tangential orbits) . 

With these assumptions, the model has 4n -I- 1 free pa- 
rameters, where n is the number of adopted components: 
namely the luminosity Ltot, scale length (r^ or h), mass-to- 
light ratio M/L and fiattening 6/a of each component plus 
the inclination angle of the galaxy. In principle, photometry 
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Table 14. Parameters from the dynamical models 



Object 




Bulf 


;e 






Disc 




Dark halo 


Bulge 4- 


Disc 




L 


Lb 


M 


M/L 


L 




M 


MIL 


M d 


M 


MIL 




[%] 


[L0,v] 


[Mq] 




[%] 


[L0,v] 


[Me] 


[To] 


[Mq] 


[M0] 


[T©] 


NGC 772 


52% 


5.8 ■ 10^0 


3.8 • IQii 


6.6 


48% 


5.4 ■ lO^o 


1.8 ■ 10" 


3.3 




5.6 ■ 10" 


5.0 


NGC 3898 


67% 


l.l-lQio 


7.1 • IQio 


6.4 


33% 


4.8 ■ 10^ 


2.0 ■ 10"' 


4.2 


9.0 ■ 10^ 


9.1 ■ 10^0 


5.8 


NGC 7782 


11% 


1.03 ■ lO^o 


8.5 ■ lO^o 


8.3 


89% 


8.6 ■ lO^" 


3.6 ■ 10" 


4.2 




4.4 ■ 10" 


4.6 



from the 2D decomposition in Table |l3| 

adopting the total corrected V — band luminosities i'yi.^ derived from Table |l| 
from the dynamical models shown in Figs. ^, [l5| and p^. 
^ at the outermost observed radius measured for the ionized-gas component. 



can be used (as explained in Sect. [l.l| ) to constrain all these 
parameters except the mass-to-light ratios and the inclina- 
tion angle j, which must be derived by kinematics. 

For each given choice of the parameters above, we com- 
pute the gravitational potential of the total mass distribu- 
tion and integrate the Jeans equations for the stellar compo- 
nent, obtaining a self-consistent model of the rotation veloc- 
ity and velocity-dispersion profiles which include the asym- 
metric drift effects. 

In order to compare the observed data with the pre- 
diction of the model, we also need to reproduce the devia- 
tions of the LOSVD profiles from a pure Gaussian shape. 
In fact, in the regions where the bulge and disc luminosi- 
ties are comparable we expect that the superposition of the 
rapid rotation of the disc with the slower rotation of the 
other components will produce a clearly non-Gaussian, and 
sometimes even 2-peaked, LOSVD even assuming that each 
individual component has a Gaussian velocity distribution. 

We parametrized these deviations in terms of the usual 
Gauss-Hermite series higher moments hz and /14, that we 
obtained by means of a first-order approximation from the 
momenta of the model velocity distribution (van der Marel 
& Franx 1993; Pignatelli & GaUetta 1999). 

Finally, we convolved the results of the model with the 
seeing and took into account the instrumental setup used in 
the different observations (slit width, pixel size). The final 
model profiles are compared with the observed stellar kine- 
matic data and the best-fit model is found with the help of 
the standard reduced x'^-analysis. Masses and mass-to-light 
ratios of different components are given in Table 

Once the fit of the stellar kinematics has been per- 
formed, and the overall potential of the galaxy is known, one 
can derive the circular velocity Vc = ii(9$/9i?) directly. By 
overlaying the Vc obtained in this way (and corrected for 
inclination) on the observed gas rotation velocity, we can 
immediately notice any deviation from purely circular mo- 
tion. The evaluation of the different convolution effects is 
crucial especially within the innermost regions, where the 
rise of the velocity curve is smoothed, and the value of the 
velocity dispersion is increased, by the seeing effects. 



5 RESULTS 
5.1 NGC 772 

5.1.1 Stellar and ionized-gas kinematics 

The stellar kinematic data extend to more than 30" (5.0 
kpc) on the receding side of NGC 772 and to about 24" (4.0 
kpc) on its approaching side (Fig. |^. In the inner 4" (0.7 
kpc), the rotation velocity of stars increases to 80 km . 
At larger radii, it rises more gently out to 120 km s"^ at 
the last measured point. At the centre the stellar velocity 
dispersion shows a maximum of about 120 km s^^; away 
from the nucleus it falls off to 75 km at 4". Outside it 
remains high, ranging between 90 km s"^ and 140 km s"^. 

The stellar velocity curve and velocity-dispersion profile 
we measured along the major axis of NGC 772 are compared 
in Fig. ^ to those obtained by Heraudeau & Simien (1998). 
The two data sets cover almost the same radial region and 
have been collected along two close position angles (P. A. = 
130° and 131°). Our line-of-sight velocities agree within the 
errors with those measured by Heraudeau & Simien (1998). 
However their velocity dispersions are higher than those we 
found. This is only due to the different technique adopted in 
measuring the stellar kinematic parameters rather than to 
real kinematic features observed at the two different position 
angles. 

The ionized-gas velocity is measured out to 45" (7.6 
kpc) on the receding side and to less than 40" (6.7 kpc) on 
the approaching side (Fig. The ionized-gas rotation ve- 
locity has a steeper gradient than the stellar velocity, reach- 
ing a value of 140 km s"'^ at |r| ~ 5" (0.8 kpc), increasing 
to 175 km s^i at |r| ~ 10" (1.7 kpc) and then flattening 
out. The gas velocity dispersion strongly peaks to about 150 
km s~i in the centre; it drops rapidly to values lower than 
50 kms-i for Irl Z 5". 



5.1.2 V—band surface photometry and bulge-disc 
decomposition 

The V—band ellipse-averaged radial profiles of surface 
brightness, ellipticity, position angle and cos 46* Fourier co- 
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Figure 2. The stellar (filled circles) and ionized-gas (open cir- 
cles) kinematics measured along the major axis (P. A. = 130°) of 
NGC 772. 



efficient of NGC 772 have been measured out to 229" (38.5 
kpc) from the centre (Fig. ^ . 

In the inner 60" (10.1 kpc) the values of ellipticity and 
position angle show a series of bumps and wiggles, oscil- 
lating between 0.15 and 0.36 and between 98° and 143°, 
respectively. These features are due to the inner spiral arms 
of the galaxy. The bump in the surface-brightness profile at 
90", and the corresponding abrupt variations in the position 
angle and fourth-order cosine Fourier coefficient can be as- 
cribed to the prominent northern arm (see Fig. ^ . At larger 
radii, the galaxy light is dominated by the contribution of 
the disc component characterized by an exponential surface- 
brightness profile and almost constant ellipticity (e ~ 0.35) 
and position angle (P. A. ~ 120°). 

In Fig. 1^ our data are compared to those obtained by 
Heraudeau & Simien (1996) in the same band. The photo- 
metric profiles of Heraudeau & Simien (1996) extend out to 
180" from the centre and their surface brightness is ^ 0.4 
mag arcsec"^ fainter than our. In the inner 50" their posi- 
tion angle and ellipticity data have smoother radial trends 
than those we measured. We noticed that the position-angle 
values included by Heraudeau & Simien (1996) in Table 5 of 
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Figure 3. NGC 772 major-axis stellar kinematics. The stellar 
kinematics derived in the present study at P. A. = 130° {filled 
circles) is shown superimposed on that obtained by Heraudeau 
& Simien (1998) at P. A. = 131° (open squares). 
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Figure 4. The V- band surface-brightness, ellipticity, position 
angle and cos A8 coefficient profiles we measured as a function 
of the semi-major axis length for NGC 772 (open circles). Error 
bars smaller than symbols arc not plotted. The V- band surface- 
brightness, ellipticity and position angle profile from Heraudeau 
& Simien (1996, filled circles) are also plotted. 



their paper differ from those they plotted in their Fig. 10. 
In our comparison we adopted the plotted ones. The dis- 
agreement between the two sets of data can be explained 
if the tabulated values result from an isophotal ellipse fit- 
ting performed on an image of NGC 772 with an incorrect 
orientation. 

In Fig. 1^ we show the result of the bulge-disc photo- 
metric decomposition of the surface-brightness distribution 
of NGC 772, which has been performed with the parametric 
two-dimensional technique discussed in Sec. 4.1. The spiral 



arms extending into the very centre of NGC 772 are clearly 
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Figure 5. The residual image of NGC 772 obtained by subtract- 
ing the model surface brightness to the observed one. The ori- 
entation of the image is the same of Fig. H. The photometric 
parameters of the model are given in Table ha. 

visible in the residual image obtained by subtracting the 
model surface brightness of the galaxy from the observed 
surface brightness. These structures resemble those detected 
by CaroUo et al. (1997, 1998) in several Sa - Sbc galaxies, 
in which the spiral pattern is visible down to the innermost 
radius accessible to the HST WFPC2 imaging. 

By fitting an exponential profile to the ellipse-averaged 
surface-brightness profile of NGC 772 measured in the B 
band, Lu (1998) derived the following disc parameters /xq = 
22.01 mag arcsec"^, h = 72'.'4 and b/a = 0.58. Our exponen- 
tial disc results larger {h = 96'.' 3) and slightly more inclined 
(6/a = 0.39) than that obtained by Lu (1998). This is due 
to the different decomposition technique we adopted more 
than to the different bandpass in which we observed. In fact 
Lu (1998) fitted his exponential disc directly to the observed 
surface-brightness profile in the radial range between 75'.' 1 
and 133'.' 1 (which he judged by eye to be disc-dominated) 
without taking into account any contribution from the con- 
spicuous bulge component. 

5.1.3 lonized-gas distribution 

In NGC 772 most of the Hii regions visible in our Hq:+[Nii] 
image (Fig. ^ lie along the two arms, extending to the north- 
ern side of the galaxy, and along the short and double-ended 
arm opposite to them. 

The ionized-gas emission is concentrated in the region of 
the galaxy (r 40") bounded by the pseudo-ring formed by 
the two tightly- wound northern spiral arms. In the Ha-f [N ii] 
emission map it is possible to disentangle an outer region ex- 
tending between 10" and 40" and characterized by smooth 
and diffuse emission, from the inner region (r ^ 10") char- 
acterized by more concentrated emission and contributing 
about 5% of the total flux. 

The ionized gas between 10" and 40" shows the kine- 
matic behaviour typical of a disc {V/a ~ 5) as inferred from 
its high rotation velocity and low velocity dispersion (Fig. 
i- 




-150 -100 -50 50 100 150 
r [arcane] 

Figure 6. The stellar-continuum image {upper panel) and the 
continuum-subtracted Ha-|-[Nll] emission image {lower panel) of 
NGC 772. 



In the inner region the kinematic behaviour of the gas is 
different. The observed gas velocity dispersion (ranging be- 
tween 70 and 150 km s~^) is comparable to that of the stellar 
component and therefore is far higher than that expected 
from thermal motion or small-scale turbulence {a ^ 50 
km s~^) typically observed in a gaseous disc. This may well 
be an indication that in this region the gas is being dynam- 
ically supported by pressure (i.e. random motions) rather 
than rotation, as it occurs in the outer part of the galaxy. 
There is, however, a wide region (4" < r < 10") in which 
the gas velocity dispersion is still only marginally greater 
than 50 km s~^, and appears to be lower than its stellar 
counterpart (see panel (b) in Fig. ^). The high velocity dis- 
persions observed near the nucleus could then only be an 
observational effect due to seeing. In order to discriminate 
between the two cases, the application of a dynamical model 
is required. 
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Figure 7. Photometric and dynamical radial profiles for NGC 772 with their respective best-fit curves obtained from our models, (a): 

V— band surface briglitncss (filled symbols) as a function of radius along tlic major axis. Wc also plot the model bulge (dotted line) 
and disc (dashed lines) surface-brightness profiles, together with their sum convolved with the seeing (solid line), (b): Ellipse-averaged 
ellipticity for the V-band data (filled circles) and the photometric model (solid line), (c): Observed stellar velocity (filled symbols) with 
its associated model (solid line), and the ionized-gas velocity (open symbols). The dashed line is the circular velocity inferred from the 
dynamical model, (d): the same as (c) but for the velocity dispersion, (e): the same as (c) but for hs coefficients of the Gauss-Hermite 
expansion of the line profile of the stars, (f ): the same as (c) but for coefficients of the Gauss-Hermite expansion of the line profile of 
the stars. In all the plots the square symbols and the triangle symbols represent data derived for the approaching NW side and for the 
receding SE side, respectively. 
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5.1.4 Dynamical modelling 

In Fig. we show for NGC 772 the comparison between 
the photometric and kinematic data and the best-fit photo- 
metric dynamical models we obtained with the techniques 
discussed in Sec. ^ The mass-to- light ratios, flattening and 
masses of bulge and disc stellar components of NGC 772 
derived with the dynamical modelling described in Sec. 4.2 
are given in Table ^ 

The agreement between the stellar kinematic data and 
the predictions of the dynamical model are fairly good at 
all radii, apart from the features observed in the LOSVD 
at |r| ~ 3". At this radius we notice a sharp decrease (from 
115 km s~^ to 80 km s~^) in the observed velocity dispersion 
and very strong asymmetries in the line shape, while the ro- 
tation velocity does not show any significant variation. The 
model is unable to reproduce these abrupt but radially con- 
fined changes in the velocity-dispersion, /13 and /14 profiles. 
Either absorbing dust or nuclear peculiar motions could be 
responsible for this feature. 

It is worth noting that the model reproduces the ob- 
served gaseous and stellar rotation velocity without the need 
of introducing a dark halo. According to the photometric 
decomposition, our kinematic data are limited to the bulge- 
dominated region of NGC 772 where the luminous mass is 
expected to dominate the galaxy dynamics. In agreement 
with these results, the dynamical model shows that asym- 
metric drift has a sizable effect on the stellar component in 
the whole observed region as it can be desumed by consider- 
ing that we derived a difference of 35 km s~^ between stellar 
rotation velocity and circular speed at the farthest measured 
point. 

Once the convolution with seeing and instrumental 
setup has been taken into account, the gas appears to rotate 
in almost circular orbits. We conclude that the gas rotation 
curve of this galaxy can not be classified as a slowly-rising 
one, and that the observed central rise in the gas velocity 
dispersion is due to seeing effects. 
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Figure 8. The stellar {filled circles) and ionized-gas {open 
circles) kinematic parameters measured along the major axis 
(P.A. = 107°) of NGC 3898. 



5.2 NGC 3898 

5.2.1 Stellar and ionized-gas kinematics 

The stellar kinematic parameters are observed out to 20" 
(1.7 kpc) and 30" (3.4 kpc) on the approaching and receding 
side of NGC 3898, respectively (Fig. N). The stellar rotation 
velocity increases linearly up to 130 km s~^ in the inner 10" 
(0.8 kpc), then it decreases to 90 km s~^ at the last ob- 
served radius. At the centre, the stellar velocity dispersion 
is about 220 km s~^, then increases up to a maximum value 
of 240 km s^^ at |r| ~ 1" (Fig. suggesting the presence 
of an unresolved kinematically decoupled core component 
(see Friedli 1996 for a discussion about this feature in the 
velcity-dispersion profile) . Off the nucleus, the stellar veloc- 
ity dispersion decreases to 160 km s~^ on the approaching 
side and 130 km s~^ on the receding side. 

The stellar velocity curve and velocity-dispersion profile 
measured along the major axis of NGC 3898 are compared 
in Fig. ^ to those obtained by Whitmore et al. (1984), Fill- 
more et al. (1986) and by Heraudeau et al. (1999). Although 
the velocity data agree within the errors, a difference in the 
central velocity gradients measured by the various groups is 



present. In particular, there is a velocity discrepancy (rang- 
ing between 100 and 150 km s~^) between our outer rota- 
tion velocities and those by Whitmore et al. (1984), who also 
measured a higher velocity dispersion on the receding side. 

The ionized-gas kinematic parameters extend to 90" 
(7.5 kpc) on each side of the nucleus (Fig. The gas rota- 
tion velocity has a steeper gradient than the stellar velocity, 
reaching 200 km s~^ at |r| ~ 10" and remaining almost con- 
stant at this value further out. The gas velocity dispersion 
is about 90 km s"'^ for |r| ^ 10" decreasing at larger radii 
(|r| ^ 30") to 20 km s^^ on both sides of the nucleus. 

The ionized-gas rotation curves we measured along the 
major axis of NGC 3898 from the spectra obtained with 
INT and MMT are compared in Fig. ^ to the Ha rotation 
velocities obtained by Rubin et al. (1985) and Fillmore et al. 
(1986). There is a good agreement between the line-of-sight 
velocities obtained in the different runs. In particular, all the 
rotation curves closely match each other, showing the same 
wiggles and bumps in the velocity for \r\ < 8". However, 
between 8" and 20" on both sides of the nucleus Rubin et al. 
(1985) measured a slower rotation than either we or Fillmore 
et al. (1986) derived. As far as the gas velocity dispersion is 
concerned, only our own data are available. They are based 
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Figure 9. NGC 3898 major-axis stellar kinematic parameters. 
The stellar velocity and velocity dispersion derived in the present 
study at P. A. = 107° {filled circles) are shown superimposed on 
those obtained by Whitmore et al. (1984, open circles), Fillmore 
et al. (1986, crosses), and Heraudeau et al. (1999, open squares) 
for the same position angle. 




17 
IB 

19 j 

o 

22 % 

23 ^ 
24 
25 
26 

25 50 75 100 125 150 25 50 75 100 125 150 
r [arcsec] 

Figure 11. The V— band surface-brightness, ellipticity, position 
angle and cos 49 coefficient profiles we measured as function of 
the semi- major axis length for NGC 3898 (open circles). Error 
bars smaller than symbols are not plotted. The band surface- 
brightness, ellipticity and position angle profiles from Heraudeau 
& Simien (1996, filled circles) and the surface-brightness profile 
by Watanabe (1983, crosses) are also shown. 
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Figure 10. NGC 3898 major-axis ionized-gas kinematic param- 
eters. The ionized-gas velocities derived in the present study at 
P. A. = 107° using the spectra obtained at INT (open circles) 
and MMT (asterisks) are shown superimposed on those obtained 
measuring the Ha emission line by Rubin et al. (1985, crosses) 
and Fillmore et al. (1986, open triangles) along the same position 
angle. 



on the line width of the [O iii] A5006.8 and Ha emission lines 
present in the MMT and INT spectra. The two data set are 
consistent within the scatter of the points derived from the 
[O III] A5006.8 line. 



5.2.2 V—band surface photometry and bulge-disc 
decomposition 

The 1/— band radial profiles of surface brightness, ellipticity, 
position angle and cos 4^ Fourier coefficient measured by 
fitting ellipses to the NGC 3898 iso pho tes extend out to 
164" (13.6 kpc) from the centre (Fig.]ll]). 

The ellipticity rises from 0.23 to 0.39 in the inner 12" 
(1.0 kpc) remaining almost constant out to 40" (3.3 kpc). 
Between 40" and 110" (9.1 kpc) it increases to a maximum 
of 0.5 at 75" (6.2 kpc) and decreases to the initial value of 
0.39. Further out it remains almost constant at 0.37 out to 
the farthest observed radius. The isophotal position angle is 
constant at 109° at all radii. The isophotes are slightly disky 
between 40" and 85" (7.1 kpc) and elliptical elsewhere. The 
surface-brightness profile is extremely regular with a few 
small bumps due to the spiral arms. 

In Fig. ^ our surface-brightness radial profile is com- 
pared to those obtained by Watanabe (1983), and Her- 
audeau & Simien (1996), who also measured for NGC 3898 
the ellipticity and position angle radial profiles. The agree- 
ment between our data and those of Heraudeau & Simien 
(1996) is good. The surface brightness measured by Watan- 
abe (1983) in the radial region between 10" and 100" is 5^ 0.5 
mag arcsec"^ fainter than ours. The surface-brightness pro- 
file of Kodaira et al. (1990) is not shown in Fig. |l^. In their 
Photometric Atlas of Northern Bright Galaxies only the pro- 
files extracted along the major and minor axes of NGC 3898 
are available and there is not an ellipse-averaged profile. This 
is also true for NGC 7782. 

The result of the two-dimensional bulge-disc paramet- 
ric decomposition of the surface-brightness distribution of 
NGC 3898 is shown in Fig. ^ The spiral pattern shows up 
in the panel of the residuals, as well as a central structure 
due to the difference in this region between the observed 
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Figure 12. The residual image of NGC 3898 obtained by sub- 
tracting the model surface brightness to the observed one. The 
orientation of the image is the same of Fig. H. The photometric 
parameters of the model are given in Table ha. 



ellipticity (increasing from 0.28 to 0.39) and the constant 
value {etuige = 0.39) derived for the bulge component. 

The surface-brightness profiles of Watanabe (1983) and 
Kodaira et al. (1990) were fitted with a r^^^ bulge and an 
exponential disc by van Driel & van Woerden (1994) and 
Baggett et al. (1998), respectively. Van Driel & van Wo- 
erden (1994) using Burstein's (1979) method obtained the 
following photometric parameters /ie = 17.1 mag arcsec"^, 
Ve = 4'.'4, /io = 20.6 mag arcsec'^ and h = SS'.'O. The 
best-fit parameters derived by Baggett et al. (1998), per- 
forming an intensity-weighted fit between 3" and 156" to 
avoid the portion of the profile most affected by seeing, are 
fXe = 18.3 mag arcsec"^, re = 6'.'2, /io = 20.8 mag arcsec"^, 
and h = 3/.' 8. The results of van Driel & van Woerden 
(1994) and Baggett et al. (1998) are consistent if we take 
into account that they used an ellipse-averaged profile and 
the profile extracted along the galaxy major axis, respec- 
tively. The bulge and disc scale parameters we obtained for 
NGC 3898 are quite different since we have taken properly 
into account the apparent axial ratios of the two compo- 
nents. 



5.2.3 lonized-gas distribution 

Ho, Filippenko & Sargent (1997) classified the nucleus of 
NGC 3898 as intermediate between an Hii nucleus and a 
LINER. This bright nucleus is also visible in our Ha+[Nii] 
emission image of NGC 3898 (Fig.|l|), mapping the ionized- 
gas distribution of the galaxy. The ionized-gas emission 
is smooth and featureless in the bulge-dominated region 
(r 20"), where about 10% of the warm gas resides. It 
shows a clumpy and fragmented distribution in the disc- 
dominated region where numerous Hii regions trace the 
NGC 3898 multiple-arm structure. The transition between 
the two regions occurs abruptly at 35" from the centre where 
the gas emission in the major-axis spectrum (Fig. ^ cuts off. 
The distribution of the H ii regions (which closely follows the 
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Figure 13. The stellar-continuum image (upper panel) and the 
continuum-subtracted Ha-|-[Nll] emission image [lower panel) of 
NGC 3898. 



continuum isophotes) indicates that the gas and stellar discs 
are coplanar. 

On the SW side of the galaxy two series of emission 
knots are aligned at a projected angular distance of 80" 
and 180" from the centre, respectively (Fig. |l^). These Hii 
regions trace the two external faint arms located far beyond 
the galaxy main body as discussed in the Carnegie Atlas of 
Galaxies (Sandage & Bedke, 1994, Panel 79). 



5.2.4 Dynamical modelling 

In Fig. |l^ we show the comparison between the photomet- 
ric and kinematic data and the best-fit dynamical model 
obtained with the technique discussed in Sec. 4.2 . The mass- 



to-light ratios of the bulge and disc stellar component are 
shown in Table |l^. 

A number of kinematic features highlight the presence 
of a hot spheroidal component. At (rj = 12" both hs and 
show the presence of a non Gaussian LOSVD which can be 
interpreted as due to the superposition of a radiply rotating 
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Figure 14. Continuum-subtracted Hq:+[Nii] emission image of 
NGC 3898 sliowing the H II regions tracing the two external faint 
arms of the galaxy. The dashed line indicates the optical size of 
the NGC 3898 corresponding to the D25 B— band diameter given 
in Tabic |^. 



disk and an almost non-rotating bulge. The stellar rotation 
velocity is everywhere less than 2/3 of the gas velocity, as a 
result of the asymmetric drift effect. Finally, for 12" < \r\ < 
25" the stellar velocity dispersion shows a plateau at a value 
of 110 km s~^ , which is the signature of a thick, dynamically 
hot component. 

The gas rotation velocity is well approximated by the 
circular velocity of our models for |r| > 8". Out to 80" from 
the centre the gas rotation curve does not show hints of an 
asymmetric drift effect and is almost flat at the constant 
value of about 200 km s~^. Kent (1988) and Moriondo et 
al. (1998b) already pointed out that for |r| ^ 8" the gas 
rotates more slowly than expected on the basis of the stellar 
kinematics and of the photometry. The comparison of the 
observed gas kinematics with our self-consistent dynamical 
model confirms this result. Moreover, we can also conclude 
that seeing convolution and finite slit width can not account 
for this effect, since these parameters were already included 
in the computation of our modeled velocity curves. We stress 
that, for this particular object, the gas rotation curve can 
not be interpreted as a sign of 'pressure-supported' gas in the 
central regions of the galaxy. First, the observed gas velocity 
dispersion is too low to account for the large differences be- 
tween the circular velocity computed from the model and the 
observed gas rotation velocity. Moreover, the gas kinemat- 
ics show strong asymmetries in the region where the 'slowly 
rising' effect takes place, with a difference of more than 80 
km s~^ in the rotation velocities of the leading and receding 
sides. Such asymmetries are not seen in the gas velocity- 
dispersion profile. It is likely that other effects are at work 
here; however, the available data do not allow us to distin- 
guish between non-axisymmetric distortion of a gaseous disc 
and possible non-equilibrium motions of the gas. 

For this galaxy only we were not able to reproduce 
the gas rotation velocities at large radii without adding a 



dark matter halo. For sake of simplicity, we adopted the 
'maximum bulge-|-disc' paradigm, in order to minimize the 
amount of dark matter needed to fit the data. A more precise 
measurement of the mass and distribution of the dark mat- 
ter would require a best-fitting procedure of the observed 
kinematic curves varying simultaneously all the parameters 
involved, but can not be applied to this object because of 
the limited extension of the gas rotation curve (« 2h; see 
Persic, Salucci & Stel 1996 for a discussion on the subject). 
We assumed a pseudo-isothermal profile for the dark halo 
with a circular velocity, given by: 



1 — — arctan 



(- 



(12) 



Due to the limitations of the available data, we were 
unable to determine both the total mass and scale radius 
of the halo. Nevertheless, the total mass within at the at 
the outermost observed radius was quite well constrained 
and turned out to be 9 • 10® Mq , with an error of less than 
10%. We also found that the scale length of the dark halo 
is so large (rn > 65"), that in the observed radial range the 
circular velocity contributed by the halo rises almost linearly 
with radius. 

Total mass estimates and mass decompositions for this 
galaxy are already present in literature (Fillmore et al. 1986; 
Kent 1998; Moriondo et al. 1998b). Remarkably, despite the 
various different hypotheses adopted, our result for the to- 
tal luminous mass agrees within 20% with all these authors: 
when scaled to our adopted distance of D = 17.1 Mpc, Kent 
(1988) gives Af,„„ = 1.1 ■ 10^^ Mq, Fillmore et al. (1986) 
Afium — 1.1 • 10^^ M0, and the most recent estimate, given 
by Moriondo et al. (1998b) is Mi„™ = 9.2 ■ 10^° Mq, almost 
identical to our value of Afium = 9.1 • lO^'^ Mq. On the con- 
trary, the mass decomposition of the different authors differs 
from each other by large factors, being strongly dependent 
on the hypothesis adopted. Kent (1988) and Moriondo et al. 
(1998b) derive bulge masses, which are respectively 15 and 3 
times lower than the value we obtained. This is due to their 
assumption that the gas rotation velocities were tracing the 
potential everywhere. As a consequence unphysical values 
for the mass-to-light ratios of both components are found, 
as the authors correctly remarked. This underestimate of the 
bulge contribution to the total mass content leads also to an 
overestimate the disc mass (by a factor of 4, if compared 
with our results) , and the mass of the dark halo (in the case 
of Moriondo et al. 1998b). 

Fillmore et al. (1986) derived the mass-to-light ratios 
from the stellar kinematics rather than from the gas; how- 
ever, their photometric decomposition is one-dimensional 
and leads to a smaller and more compact bulge in combina- 
tion with a larger disc than the ones we derived. Therefore 
their results are not easily comparable with ours. On the 
other hand, our two-dimensional photometric decomposition 
give a bulge-to-disc ratio in good agreement with both Kent 
(1988) and Moriondo et al. (1998b), who did not use for 
their decomposition the ellipse-averaged surface-brightness 
profiles as done by Fillmore et al. (1986). 



5.3 NGC 7782 
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Figure 15. Same as Fig. |^ but for NGC 3898. The squares and triangles represent data measured on the receding SE and approaching 
NW sides, respectively. In panel (c) the model circular velocity is plotted with (dot-dashed line) and without (dashed line) the contribution 
of a dark matter halo with the structural parameters given in Table |l4l 
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Figure 16. The stellar (filled circles) and ionizcd-gas (open 
circles) kinematic parameters measured along the major axis 
(P. A. = 175°) of NGC 7782. 



5.3.1 Stellar and ionized-gas kinematics 

The stellar velocity curve is observed out to 40" (14.6 kpc) 
on both sides of NGC 7782 (Fig. The stellar rotation 
velocity increases linearly up to 70 km in the inner 1" 
(0.4 kpc). It does not change for 1" <, \r\ <, 3" (1.1 kpc), 
and further out it shows a shallower gradient rising to 240 
km s~^ at 14" (5.1 kpc). Outwards the stellar rotation ve- 
locity remains constant or possibly even rises. At the centre 
the stellar velocity dispersion peaks to about 190 km s~^, 
then off the nucleus it decreases gradually to less than 50 
km s~^ at the last outermost radii. 

The ionized-gas velocity curve extends to less than 30" 
(10.9 kpc) on the approaching side and to less than 45" 
(16.4 kpc) on the receding side of the galaxy (Fig. |l^. The 
gas rotation velocity has a steeper gradient than the stellar 
curve, rising to 240 km s~^ at 3"; then it remains almost 
constant, becoming similar to the stellar velocity for |r| > 
14". The gas velocity dispersion is centrally peaked at about 
180 km s~^, dropping to values lower than 100 km s~^ for 
radii larger than 3" on both sides of the galaxy. Further 
out it oscillates around values of about 40 km s~^. The gas 
velocity dispersion is not greater than the stellar dispersion 
over the whole observed radial range. 




Figure 17. The V- band surface-brightness, ellipticity, position 
angle and cos 46 coefficient profiles as function of the semi-major 
axis length for NGC 7782. Error bars smaller than symbols are 
not plotted. 



5.3.2 V —hand surface photometry and bulge-disc 
decomposition 

The V— band radial profiles of surface brightness, ellipticity, 
position angle and cos 49 Fourier coefficient of NGC 7782 
are measured out to 83" (30.3 kpc) from the centre (Fig. 

The ellipticity increases from 0.18 to 0.37 in the first 10" 
(3.6 kpc), dropping to 0.22 further out. It then increases 
again to 0.48 and flattens around this value for r > 30". 
The position angle is constant at 175° at almost all radii, 
except for two abrupt changes at 6" and 18" by less than 20° . 
Isophotes are boxy-shaped around 6" and between 18" and 
36", and almost elliptical elsewhere. The surface-brightness 
profile is characterized by a plateau at 21.6 mag arcsec"'^ 
between 23" and 30". All these features can be explained as 
due to the spiral pattern (see Fig. |l8|) . 

The band surface-brightness radial profile of 
Koidaira et al. (1990) has been extracted along the galaxy 
major axis but it is not directly comparable with the ellipse- 
averaged one we have derived and shown in Fig. 

The two-dimensional bulge-disc parametric decompo- 
sition of NGC 7782 was performed adopting an exponen- 
tial bulge since this yielded lower residuals. The map of the 
residuals obtained as a difference between the observed and 
the model surface brightness is plotted in Fig. This fig- 
ure shows the asymmetric pattern of the inner stellar spiral 
arms. 



5.3.3 Ionized-gas distribution 

The ionized-gas distribution of NGC 7782 is shown in Fig. 
p^ . The map of the Hq-|-[Nii] emission is characterized by 
the presence of a small and bright nuclear region aligned 
with the galaxy major axis, and two tightly wound spiral 
arms which can be followed for 360° forming a sort of double 
ring-like structure. The gaseous arms are more symmetric 
than the stellar arms and do not extend to the inner region of 
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Figure 18. The residual image of NGC 7782 obtained by sub- 
tracting the model surface brightness to the observed one. The 



orientation of the image is the same of Fig. 
parameters of the model are given in Table 



16| The photometric 



the galaxy, starting outside the constant surface-brightness 
region at 30" from the centre. 

As in the case of NGC 772, the bright Hq-|-[N ii] nucleus 
of NGC 7782 is characterized by a high ionized-gas velocity 
dispersion {ogas > 100 km for |r| ^ 3") as shown in Fig. 
p^ , suggesting the possible presence of pressure-supported 
gas in the bulge region. In spite of its apparent clumpy dis- 
tribution (Fig. ^9|) a diffuse component of ionized gas is also 
present in the disc of NGC 7782, as indicated by the fact that 
we are able to measure continuously the ionized-gas velocity 
parameter without gaps along the complete major axis out 
to 30" and 45" on the two sides of the nucleus (Fig. [l6| ) 
respectively. The distribution of the Hii regions, following 
the continuum isophotes, indicates that the gas and stellar 
discs are coplanar. 
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Figure 19. The stellar-continuum image (upper panel) and the 
continuum-subtracted Ha-|-[Nll] emission image (lower panel) of 
NGC 7782. 



5.S.4 Dynamical modelling 

The comparison between the photometric and kinematic 
data and the best-fit photometric dynamical models we ob- 
tained for NGC 7782 is shown in Fig. ^ The mass-to-hght 
ratios, the flattening and the masses of the bulge and disc 
stellar components derived using the dynamical modelling 
technique described in Sec. 4.2 are given in Table PL4l 



From both the kinematic and the photometric data 
shown in Fig. ^ we can clearly distinguish the bulge and 
disc-dominated regions. 

An exponential and almost spherical bulge is the main 
photometric and dynamical component within 4" from the 
centre. In this region the stellar velocity dispersion shows a 
plateau at 180 km s~^ and in agreement with the asymmet- 
ric drift effect calculated by the model the stellar rotation 
velocity is below the value of the circular velocity. In the 
transition region bewteen 2" ^ |r| 7" the gas rotation 
curve shows small deviations (IV^as — Vdrc] ^ 50 km s~^) 
from pure circular rotation. Any evidence for possible non- 
circular motions is marginal and the increase of the gas ve- 



locity dispersion towards the nucleus can be better explained 
as due to a seeing effect. 

For \r\ > 15" the surface brightness and ellipticity pro- 
files of NGC 7782 are dominated by the light contribution of 
the disc component. Therefore in the outer region asymmet- 
ric drift effect is negligible and both gas and stars appear 
to rotate close to the circular velocity. No dark matter halo 
is required to reproduced the observed velocity curve out to 
the last observed radius. 

Since the /13 and /14 values are poorly constrained by 
observations, it is not possible to draw any conclusion about 
anisotropics in the stellar velocity dispersion. 



6 DISCUSSION AND CONCLUSIONS 

Non-circular velocities of the ionized gas in the bulge- 
dominated region of early-to-intermediate type disc galax- 
ies have been detected by different authors (Fillmore et al. 
1986; Kent 1988; Kormendy & Westpfahl 1989). The discov- 
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Figure 20. Same as Fig. |j but for NGC 7782. The squares and triangles represent data measured on tlie receding NW and approaching 
SE sides, respectively. 
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ery that in a large fraction of SO's the gas velocity dispersion 
remains as high as the stellar velocity dispersion over an ex- 
tended radial range (Bertola et al. 1995) corroborated the 
early suggestions that at small radii gas could be dynami- 
cally hot rather than settled in a cold disc. In this framework, 
the dynamical modelling by Cinzano et al. (1999) showed 
that in the SO galaxy NGC 4036 the difference between the 
circular velocity curve inferred from stellar kinematics and 
the rotation curve of the gas is substantially accounted for 
its high velocity dispersion. 

It is usually customary, in measuring mass distribution, 
to assume that the ionized gas is moving in circular orbits. It 
is therefore crucial to understand if and when the warm gas 
in galaxies has a rotation curve which rises more slowly than 
the gravitational equilibrium circular velocity. Ignoring this 
effect can lead to a significant underestimate of the masses 
of the bulge, disc and dark halo (e.g. Kent 1988 for a dis- 
cussion). Up to now such a direct comparison between the 
circular velocity inferred from dynamical modelling based 
on stellar kinematics and surface photometry, and the ob- 
served gas rotation velocity has been done for only a small 
set of lenticulars and early-to-intermediate spirals (see Ta- 
ble |l|). The three new cases discussed in this paper repre- 
sent therefore a useful step in understanding how common 
the phenomenon of slowly-rising rotation curves is in bulge- 
dominated galaxies. 

For NGC 772, NGC 3898 and NGC 7782 we present a 
self-consistent Jeans model for the stellar kinematics, adopt- 
ing the light distribution of bulge and disc derived by means 
of a two-dimensional parametric photometric decomposition 
of the observed 1/— band surface brightness. This allows us 
not only to investigate the presence of non-circular gas mo- 
tions, but also to derive the mass distribution of luminous 
and dark matter in these objects. 

In NGC 772 the observed velocity dispersion of the ion- 
ized gas is comparable to the stellar velocity dispersion and 
is far higher than expected from the thermal motions or 
small-scale turbulence, typically observed in a gaseous disc. 
However, dynamical modelling showed that the rotation of 
the ionized-gas traces the circular velocity and the central 
rise observed in its velocity dispersion is due to seeing effects, 
which serves as a warning against over interpretation. 

On the contrary, the gas rotation curve of NGC 3898 
can be classified as 'slowly-rising', since for the inner 8" ion- 
ized gas is rotating more slowly than the circular velocity 
predicted by dynamical modelling. NGC 3898 is the only 
galaxy of our sample for which previous determinations of 
the masses of the bulge and disc components are available in 
the literature (Kent 1988; Fillmore et al. 1986; Moriondo et 
al. 1998b). Comparison shows that NGC 3898 is a clearcut 
example that mass decomposition based only on emission- 
line rotation curves can be unreliable, at least for bulge- 
dominated galaxies. This error was made and discussed by 
Kent (1988) and Moriondo et al. (1998b), who assumed the 
gas to be rotating at circular velocity at all radii and found 
extraordinarily low values for the mass-to-light ratio of the 
spheroidal component, as if NGC 3898 was an almost 'bulge- 
less' galaxy in spite of its overall morphology (see Table 
|l| and Figs. |, and |l^) and bulge-to-disc luminosity ratio 
[B/D = 2). The asymmetry of the gas rotation curve of 
NGC 3898 (with a maximum AV between the two sides of 
about 80 km s~^ at 6" from the centre), and the strong dif- 



ference between the gas {ugas < 90 km s~^) and the stellar 
velocity dispersion {astars > 120 km s~^ with a central peak 
to 210 km s~^) in the innermost 5" (where we observe al- 
most the same velocity gradient for both gas and stars) are 
an indication that along with random motions other phe- 
nomena could also be contributing to the slowly rise of the 
gas velocity. This is the case of other galaxies with slowly- 
rising rotation curves, such as the Sa NGC 2775 (Corsini 
et al. 1999) and the SO NGC 4036 (Cinzano et al. 1999). 
The asymmetry of the inner parts of the gas rotation curve 
of NGC 2775 has been explained as due to a gaseous com- 
ponent which is not rotating in the galaxy plane, and the 
possible presence of drag forces between the ionized gas and 
the hot component of the interstellar medium of NGC 4036 
has been suggested (for a discussion of the phenomenon see 
Cinzano et al. 1999). 

Finally, in NGC 7782 the gas rotation curve shows small 
deviations from pure circular rotation only in the transition 
region between bulge and disc. Any evidence for possible 
non-circular motions is marginal and the sharp increase of 
the gas velocity dispersion towards the nucleus can be better 
explained as due to a seeing effect as in the case of NGC 772. 

As far as the presence of dark matter is concerned, we 
infer that the mass is essentially traced by light in NGC 
772 and NGC 7782, where gas rotation velocities were ob- 
served out to 0.2i?25 and 0.6i?25, respectively. For NGC 
3898 the combined stellar and gaseous rotation data (which 
extend out to 0.2_R25 and 0.7-R25 respectively) require the 
presence of a massive dark halo. This result is based on the 
idea that gas kinematics at large radii is representative of 
a dynamically cold disc supported by rotation, as suggested 
by the high gas rotation velocity and low velocity dispersion 
(V/cr ~ 10) for r > 30" (corresponding to about 0.2 R25). 
This result also qualitatively agrees with the general dark 
matter scenario, with the less massive galaxies being the 
more dark-matter dominated (Salucci & Persic 1999), al- 
though these data do not permit a careful measurement of 
the dark halo properties. 

Kent (1988) found that, out of the 14 Sa galaxies he 
modeled, six show gas rotation curves that rise too slowly for 
a constant mass-to-light ratio implying they are respectively 
too low for the bulges and too high for the disc components. 
According to Table |l^, the gas motion is not circular in the 
inner regions of 6 of the 9 early-to-intermediate disc galaxies, 
for which dynamical modelling allows the direct comparison 
between the gaseous and circular speeds (among these only 
NGC 3898 belongs also to Kent's sample). 

Without pretending to draw any statistical conclusion 
based on so small number of objects, it nevertheless seems 
probable that non-circular gas velocities could be a com- 
mon feature in the central kiloparsec of disc galaxies rang- 
ing from SO's to Sb's. However, according to the available 
data, it is not possible to derive a straightforward correla- 
tion with the galaxy morphological type, total luminosity, 
bulge-to-disc mass ratio or bulge size. It seems that even 
if the gas may be supported by non-circular random mo- 
tions, this type of dynamical-pressure support may not be 
the only effect for the 'slowly-rising' gas rotation curves. In 
fact, the central gas velocity dispersion is too low to account 
for the difference between gas and circular speeds in two 
of the three galaxies (namely NGC 2775, and NGC 3898) 
shown in Table O for which the gas velocity falls short of 
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Table 15. Early-to-intermediate type disc galaxies with available comparison between gas rotation velocity and circular speed 



Object 


Morphological Type 


D 




^hulge Iodise 








rv 


ra 


Ref. 




[RSA] 


[RC3] 


[MpcJ 


[mag] 




[kpcj 


[kpcJ 




[kpcJ 


[kpcJ 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


NGC 772 


Sb(rs)I 


.SAS3.. 


34.7 


-21.61 


2.11 


11.0 


6.2 


no 




0.5 


1> V 


NGC 2179 


Sa 


.SASO.. 


33.7 


-19.81 


2.80 


2.0 


5.7 


no 




0.2 


2, R 


NGC 2775 


Sa(r) 


.SAR2.. 


13.9 


-19.91 


1.39 


3.6 


4.4 


yes 


2.0 




2, r 


NGC 3898 


Sal 


.SAS2.. 


17.1 


-19.84 


3.55 


1.6 


2.3 


yes 


0.7 




1, V 


NGC 4036 


S03(8)/Sa 


.L..-.. 


20.1 


-20.03 


2.04 


1.2 


0.7 


yes 


1.0 


0.4 


3, y 


NGC 4450 


Sab pec 


.SAS2. 


17.0 


-20.40 


2.69 


1.7 


1.4 


yes 


> 5.0 


? 


A,B 


NGC 4569 


Sab(s)I-II 


.SXT2.. 


17.0 


-21.36 


0.08 


0.3 


0.4 


yes 


2.6 


? 


A,B 


NGC 5055 


Sbc(s)II-III 


.SAT4.. 


7.5 


-20.35 


0.12 


0.4 


0.3 


yes 


0.5 


? 


4,B 


NGC 7782 


Sb(s)I-II 


.SAS3.. 


74.8 


-21.95 


0.24 


1.1 


1.5 


no 




0.4 


1, V 



NOTE — col. 2: morphological type from RSA; col. 3: morphological type from RC3; col. 4: distance derived as Vo/Hq with 
Vb the velocity relative to the ccntroid of the Local Group obtained from the heliocentric systemic velocity as in RSA and 
Ho = 75 km s-i Mpc-l. The hcUoccntric velocity of NGC 5055 has been taken from RC3. For NGC 4450 and NGC 4569 
which belong to the Virgo Cluster the distance has been taken following Frecdman ct al. (1994); col. 5: absolute corrected 
B magnitude from in RC3; col. 6: bulge-to-disc mass ratio; col. 7: effective radius of the r^/'* bulge. For NGC 7782 
the bulge has an exponential surface-brightness profile; col. 8: radius at which Hbulge = lJ-disc\ col. 9: yes = the ionizcd-gas 
component shows a slowly-rising rotation curve (i.e. Vgas < Vdrc at small radii), no = the gas rotation traces circular speed 
at all radii {Vgas = Vdrc at all radii); col. 10: radial range in which Vgas < Vdrc', col. 11: radial range in which Cgas — Cstars, 
? = the radial profile of the ionized-gas velocity dispersion is not available; col. 12: references and bands of the data [1 = 
this paper, 2 = Corsini et al. (1999), 3 = Cinzano et al. (1999), 4 = Fillmore et al. (1986)]. 



the circular speed and the velocity-dispersion profile of the 
gaseous component is available. Dynamical modelling based 
both on stellar and gaseous kinematics of a larger number of 
bulge-dominated galaxies is needed to understand the links 
between slowly-rotating gas and bulge properties, as a first 
step to put constraints on its physical nature (e.g. Mathews 
1990) and origin (e.g. Bertola et al. 1995). 
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